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Abstract 

The main substance of the paper concerns the growth rate and the 
classification (ergodicity, transience) of a family of random trees. In 
the basic model, new edges appear according to a Poisson process of 
parameter A and leaves can be deleted at a rate p. The main results 
lay the stress on the famous number e. A complete classification of the 
process is given in terms of the intensity factor p = X/fj,: it is ergodic 
if p < e -1 , and transient if p > e _1 . There is a phase transition 
phenomenon: the usual region of null recurrence (in the parameter 
space) here does not exist. This fact is rare for countable Markov 
chains with exponentially distributed jumps. Some basic stationary 
laws are computed, e.g. the number of vertices and the height. Various 
bounds, limit laws and ergodic-like theorems are obtained, both for 
the transient and ergodic regimes. In particular, when the system is 
transient, the height of the tree grows linearly as the time t — > oo, at 
a rate which is explicitly computed. Some of the results are extended 
to the so-called multiclass model. 



Key words Random trees, ergodicity, transience, nonlinear differential 
equations, phase transition 



*Guy .Fayolle@inria.fr - INRIA - Domaine de Voluceau, Rocquencourt BP 105, 
78153 Le Chesnay Cedex, France 

^krikunOlbss . math . msu . su - Laboratory of Large Random Systems - Faculty of Math- 
ematics and Mechanics, Moscow State University, 119899, Moscow, Russia 

* Jean-Marc.Lasgouttes@inria.fr - INRIA - Domaine de Voluceau, Rocquencourt 
BP 105, 78153 Le Chesnay Cedex, France. 

§J.-M. Lasgouttes worked partly on the present study while spending a sabbatical at 
EURANDOM in Eindhoven. 



1 



1 Introduction and model description 



So far, very few results seem to exist for random trees as soon as insertions 
and deletions are simultaneously permitted (see e.g. |14|). We shall study 
one of the simplest models in this class, which offers both interesting and non 
trivial properties. Broadly speaking, one might think of a vertex as being a 
node of a network (e.g. the Internet) or of some general data structure. This 
paper is a self-contained continuation of the report [5]. 

Let G = {G(t),t > 0} be a continuous time Markov chain with state space 
the set of finite directed trees rooted at some fixed vertex vq- 

Throughout the study, the distance between two vertices is the number of 
edges in the path joining them, and the height h(v) of a vertex v is the 
distance from the root. The set of vertices having the same height k form 
the k-th level of the tree, the root vo being at level 0. Hence the height of 
G is a stochastic process {Hc(t),t > 0}, where 

Hc(t) = max h(v). 

veG(t) 

Na(t) will stand for the volume of G(t) (i.e. its total number of vertices). 

Wherever the meaning is clear from the context, we shall omit the subscript 
G and simply write H or N. The indegree of a vertex v is the number of 
edges starting at v and a vertex with indegree is a leaf. Finally, we will also 
need the classical notion of subtree with root v, which goes without saying. 

At time t = 0, G(0) consists of the single vertex Vq. Then at time t > 0, the 
transitions on G are of two types: 

• Adjunction. At each vertex v, a new edge having its origin at v 
can be appended to the tree at the epochs of a Poisson process with 
parameter A > 0. In this case, the indegree of v is increased by one 
and the new edge produces a new leaf. 

• Deletion. From its birth, a leaf (but the root) can be deleted at a 
rate \i. In other words, a vertex as long as it has no descendant has 
an exponentially distributed lifetime with parameter /U > 0. 



2 



1.1 Organization of the paper, results and related studies 

Section [2] is devoted to the birth and death model described above, with 
A, fi > 0. An exact and complete classification of G is given. Indeed, nec- 
essary and sufficient conditions are derived for the process to be ergodic 
(fji > Ae) or transient (fi < Ae). A phase transition phenomenon is en- 
lightened, which corresponds precisely to the absence of a null recurrence 
region. 

When the system is ergodic, the stationary distributions of the volume and 
of the height of the tree are computed in section [3l 

Section 0] deals with limit laws and scalings for Hc{t) and Na(t) in the 
transient case. The main outcome is a kind of ergodic theorem for Hc(t), 
valid for any [i > 0. It allows, in the particular case [i = (pure birth- 
process), to rediscover the magic growth rate Ae, originally derived e.g. in 

mm- 

Finally, section [5] proposes an extension to a multiclass model, in which the 
parameters of the process depend possibly on the class, the key result being 
a qualitative theorem for ergodicity. 

Recently, the authors were made aware of a model studied in [16] [13]. The 
setting considered there is a contact process on a <i-ary ordered tree, also 
known as a Catalan tree, where d > 2 is an arbitrary finite integer. The main 
difference with our model resides in the fact that each empty descendant of 
an occupied vertex can become occupied at a rate (3 > 0. The classification 
of the process was obtained for d = 2 in [16] . and |13| extends the result 
to any finite d, but nothing was said for d = oo. It turns out that most of 
the points presented in our study (ergodicity, zero-one laws, etc.) cannot 
be obtained by simply letting d — > oo in |16[ I13j . Likewise, reversibility 
arguments (which should theoretically lead to explicit invariant measures) 
used in the latter papers do not seem to be effective when d is infinite (see 
section [3]) . 

2 The birth and death case: A > 0, fi > 

The random tree G evolves according to the rules given in the introduction, 
the first important question being to find exact conditions for this process 
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to be recurrent or transient. Main results in this respect are stated in theo- 
rem [2TTJ 

For convenience, we define the lifetime t v of an arbitrary vertex v, which 
measures the length of the time interval between the birth and the death of 
v (for consistency r v = oo if v is never erased) . 

Lemma 2.1. All vertices, but the root, have the same lifetime distribution 
p(t), which satisfies the following system (S) 

P(t) = fMexp{-\j\l-p(x))dx}, (2.1) 

m = ^Ir + fiM-vWPM, (2-2) 
with the initial condition p(0) = 0. 

Proof. Let v be a particular vertex of G{t) and consider the related random 
subtree with root v. Its evolution does not depend on anything below v, as 
long as v exists. Therefore all these subtrees are identically distributed and, 
accordingly, their vertices have the same lifetime distribution. 

To capture more precisely the evolution of the process, we associate with 
each vertex v with age t its number X v (t) of direct descendants (i.e. who are 
located at a distance 1 from v). 

At rate A, a vertex v produces descendants whose lifetimes are independent, 
with the common distribution pit). As soon as X v (t) = 0, v can die at rate 
//, in which case the process of production stops. 

It is actually useful to extend X v (t) for all t > by deciding that, instead of 
deleting v, a fj,- event occurs without stopping the production of descendants. 
With this convention, the number of descendants of the root vertex vq evolves 
as X Vo (t), for all t > 0. Let t v denote the random epoch of the first //-event, 
which is distributed according to p(t). 

Clearly the process X v is regenerative with respect to the //-events. Thus 
the random variables X v (t) and X v (t v + 1) have the same distribution. 

For any fixed t, we write down a sum of conditional probabilities, expressing 
the fact that v had exactly k descendants, who all have died in [0, t], their 
birth-times being independent and uniformly spread over [0,t]. This yields 
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at once equation (|2.ip . since 



k=0 ' K { J 

t 

= expj - Ay (1 -p(s))dx}. (2.3) 



By means of a regenerative argument, it is also possible to rewrite the above 
probability in another way, starting from the decomposition 

P{X v (t) = 0} = P{X v (t) = 0,r v >t} + P{X v (t) =0,t v < t}. (2.4) 

In fact, we have the trite relations 

r ^B = ti P{x v (t) = o,r v >t}, 

at 

P{X v (t) =0,T v <t} = P{X v (t - T V ) = 0,T v <t}, 

P{t v € (y,y + dy)} = dp(y), 
which yield in particular, 

P{X v (t) = 0,T v <t}= [ P{X v (t-y) = 0}dp(y). 

Jo 

Hence, putting f3{t) = fiP{X v (t) = 0}, one sees that (|2.4|) corresponds term 
by term to (|2.2p . The proof of the lemma is concluded. ■ 

It is convenient to introduce now r = t Vq , which is the random variable 
representing the epoch of the first //-event for the root of the tree. Nonethe- 
less in the sequel, especially in sections [3] and HI r will also often refer to 
the lifetime of an arbitrary generic vertex v, owing to the fact that all these 
quantities have the same distributions. 

We are ready to state the main result of this section. 
Theorem 2.1. 

(A) The Markov chain G is ergodic if, and only if, 

def A 1 . . 

p = -<-. (2.5) 
/x e 
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(B) When the system is ergodic, the mean lifetime m = Er is given by 

r 

where r < 1 denotes the smallest root of the equation 

re~ r = p (2.6) 

and represents the mean number of descendants of an arbitrary vertex 
at steady state. 

(C) When p > — , then the system is transient. In this case, 

lim p(t) =£<1. 

t— >oo 

As a rule, x being the positive root of xe x = p -1 , we have for any p 

x < i < minfl, — \ and lim pi = 1. 

V p) p— >oo 

The proof of the theorem is spread over the next two subsections. 
2.1 Ergodicity 

Relying on standard theory of Markov chains with countable state space 
(see [6l vol. I]), we claim the system ergodic if, and only if, m < oo. As 
a matter of fact, the //-events are regeneration points for the process X(t), 
which represents exactly the number of descendants of the root vq. Hence 
when Er < oo (i.e. /3(oo) > 0), the event {X(t) = 0} has a positive prob- 
ability so that G is ergodic. Conversely, if Et = oo then X(t) is transient 
and so is G. 

For an arbitrary positive function /, denote by /* its ordinary Laplace trans- 
form 

/•oo 

/*( s )H f / e-"f(t)dt, 5R(s) > 0. 
Jo 

Later on we will also need the associated inversion formula (see e.g. [7]) 

-i /•er+ioo 

f{t) = ^ e st P{s)ds, R((t)>0. (2.7) 

^ i7r J (T — lOO 
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To show the necessity of condition (|2.5|) . we suppose G is ergodic. In this 
case, by (|2.ip . the quantity lim f3(t) does exist and 

t— >oo 

so that we can apply a limiting relation of Abelian type (see e.g. [7]). Hence 
equations (|2.ip and (|2.2p — the latter belonging to the Volterra class — yield 
respectively 

lim Pit) = lim s(3*(s) = lim P K - = — , 
t t->co s^o s-s>o 1 — sp*(s) m ^ 8) 

lim P(t) = fie~ Xm , 

V t— >oo 

whence the equality 

p = Ame- Am . 

As the function xe~ x reaches its maximum e _1 at x = 1, we conclude that 
necessarily p < e . 

In order to prove the sufficiency of (|2.5p , we have to get a deeper insight into 
system (S). There will be done along two main steps. 



(a) Although (S) reduces to a second order nonlinear integro-differential 
equation, this does not help much. What is more useful is that all deriva- 
tives p( n \0), f3^ n \0), taken at the the origin in the complex t- plane, can be 
recursively computed for all n. This can be checked at once, rewriting (|2.ip 
in the differential form 

^ + A(l-pP(i) = 0. (2.9) 

Noticing the derivatives p( n \0) and fi( n \0) have alternate signs when n 
varies, it is direct to verify that f3 and p are analytic functions around the 
origin, and that their respective power series have a non-zero radius of con- 
vergence. The first singularities of p and f3 are on the negative real axis, 
but not easy to locate precisely. Thus (S) has a solution, which is unique, 
remarking also that uniqueness is a mere consequence of the Lipschitz char- 
acter of dp{t) jdt with respect to f3 in the Volterra integral equation (|2.2p (see 
e.g. [2]). En passant, it is worth noting that the solution in the whole com- 
plex plane — which is not really needed for our purpose — could be obtained 
by analytic continuation directly on system (S). 
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(b) When (|2.5p holds, the next stage consists in exhibiting a non-defective 
probabilistic solution p(t) [necessarily unique by step (a)], with a finite mean 
m < oo. This is more intricate and will be achieved by constructing a 
converging iterative scheme 

Consider the system 
>o(i) 
Pk(t) 

< 

Pfc(O) 

The second equation in fj2. lOj) is equivalent to 

allowing to derive pk from /3k by means of ([2.7p (see also [6] for various 
inversion formulas in the real plane). Then the computational algorithm 
becomes simple: 

1. p (t) = 1 - e~i*. 

2. Compute f} x (t) = /uexp[-p(l - e - ' 4 *)]. 

3. Compute pi(t), then foit) , P2(t) , etc. 

In the scheme (|2.10p . the initial condition /3o(t) = fj, is tantamount to take 
an implicit fictitious function, say p~i, satisfying p-i(t) = 1, Vi > 0. 

At each step, the successive p^s are non-defective probability distributions, 
with finite means denoted by m&. The scheme (I2.10j) enjoys two nice prop- 
erties. 

(i) It is monotone decreasing. Suppose Pk(t) < pk-i(t), which is in particular 
true for k = 1. In the third equation of (|2,10p . is simply the 

probability of being empty for an m/g/oo queue with arrival rate A and 
service time distribution function p^. It is therefore possible, by a coupling 
argument, to build two m/g/oo queues corresponding to /3}-(t) and /3k+i(t) 
such that the //-event pertaining to level k + 1 will always occur later than 
the one for level k. Thus pk+i(t) < Pkif). 



fj,, t>0, 

^ + J%k(t-y)dp k (y), 



/iexp 
0, Vfc > . 



{-a/(i-p* 



(2.10) 
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So, the positive sequences {pk(t), (3k(t),k > 0} are uniformly bounded and 
non-increasing for each fixed t. Consequently, 

p(t) = lim\ p k (t) and p{t) = lim\ f3 k (t) 
form the unique solutions of (S). 

(ii) Letting r k = Am/% and combining the two main equations of (|2.10p . we 

get 

r k+i = P& rk > Vfc > 0, with ro = p. 

When p < e -1 , the r^'s form an increasing sequence of positive real numbers, 
with a finite positive limit r satisfying equation (|2.6p . Since 1— Pk(t) is also an 
increasing sequence of positive functions, the theorem of Beppo Levi ensures 
the equality 

roc 

(l-p(t))dt= lim / (1 -p k (t))dt = lim m k = -. (2.12) 

fe— >oo Jo k—^oo A 

It is worth to point out that (|2.10p is equivalent to the construction of a 
sequence of trees {G k , k > 0}, such that, for any finite k, G k is ergodic and 
has a height not greater than k. 

This completes the proof of points (A) and (B) of the theorem. ■ 



2.2 Transience 

It turns out that the classification of the process for p > e _1 can be obtained 
from analytic arguments. 

Recalling that i = lim p(t), we define 



t— too 

dof , / / „ / \ \ t i • ^ / , \ clef 



e(t) = X (£ — p(x))dx, lim | e(t) = e, (2.13) 



o 



and it will be convenient to write e'(t) 



t—t-oo 
def de(t) 



The quantity e exists, but a priori is not necessarily finite. It represents the 
mean number of descendants of an arbitrary vertex v, conditioning on the 
fact that v is the root of an almost surely finite tree. In fact we will show 
that if p > e~ 1 then £ < 1, in which case the system is transient. Thus 
there is no null-recurrence region in the parameter space. However, e given 
by (|2.13p will appear to be finite for all p > 0. 
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Unless otherwise stated, s in this subsection will stand for a positive real 
variable. Then by a direct computation we get 



POO 

/ exp[-(e(*) + (A(l -*) + «)*)] eft, 
Jo 



together with the functional equation 

s 2 e*(s) £ + (£-l)p*(s) 



(2.14) 



(i) Suppose now for a while £ = 1. Then (|2.14p reduces to 

s 2 e*(s) [l + f3*(s)] = A, (2.15) 

and the idea is to show that, for p > e _1 , (|2,15p has no admissible solution, 
i.e. a solution such that \im\e'(t) = 0. 

t— >oo 

By well known theorems for Laplace transforms (see e.g. [7]), we have the 
relations 

lim se*(s) = e, lim s 2 e*(s) = 0, lim sf3*(s) = pexp(-e). 
s— >o s— >o s— s>0 

When the system is ergodic, e < oo and the above limit equations give at 
once 

lim s 2 e*(s)l3*(s) = ^eexp(-e). (2.16) 

In the case e = oo, the question is more difficult and (|2.16p does not hold 
without additional conditions on e(t), as for instance slow variation (see 
tauberian theorems in |6J). The only cheap by-product of (|2.15p is the exis- 
tence of the decreasing limit 

lim f s 2 e*(s)/3*(s) = A - lim\ s 2 e*(s). 

To get deeper insight into f|2. 15[) . we remark that the quantity se*(s)f3*(s) 
can be viewed as the Laplace transform of a convolution measure with density 

t 

e (z) exp[— e(t — z)]dz, 

o 

so that (I2.15P is equivalent to the integro-differential equation 

1 /"* 

p = - £ '(t) + / e'(z) exp[-e(t - z)]dz. (2.17) 
A* Jo 
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It is worth remembering that we are searching for solutions of (|2.17p in the 
class of functions e(t) which have positive monotone decreasing derivatives 
[this property can in fact be established by taking derivatives of higher order 
in (|2,17p ]. and satisfy the intial condition e(0) = 0. Hence, as t — > oo, the 
limit of the integral in (|2.17|) exists and 



lim 

t— >oo 



1 z'(t) 
e'(z) exp[-e(t - z)]dz = p - lim\ — — > 0. 



t— >oo p 



Choose T, < t < T. Then the decomposition 



1 -e'(T)+ f e'(z)exp[-e(T - z)]dz + f e'(z) exp[-e(T - z)]dz 
1 Jo Jt 

yields by monotonicity 

i r T 

p < -e'(T) +s(t)exp[-e(T- t)] + / e'(z) exp[-e(T - z)\dz. (2.18) 
M Jt 

Putting T = 2t in (|2.18p and using e'(2t) < e'(t), we obtain the main 
inequality 



p - e(t) exp[-e(t)] < e'{t) 



1 /■* 

— h / exp[— e(z)]dz 
V Jo 



(2.19) 



On the other hand, (|2. lTj) shows immediately that the right-hand side mem- 
ber of (I2.19P is bounded by p. Finally, any solution of (|2.17p must satisfy 



p - e(t) exp[-e(t)] < e'(t) 
Assume also now p > e~ l . Then 



i r* 

— hi exp[— e(z)]dz 
/" Jo 



< 



(2.20) 



sup[p — e(t) exp[— e(t)]] = ^ > 0, 



so that (I2.20P implies 



e < e'{t) 



1 /■* 

— h / exp[— e(z)]dz 
V Jo 



(2.21) 



Let a G [£, p] be a real parameter and consider the differential equation 



f(t) 



1 f 

- + / exp[-/(z)]dz 
M Jo 



a. 



(2.22) 
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The next step is to show that all solutions of (|2.21|) are located in the strip 
delimited by the maximal and minimal solutions of (|2.22p . Setting g(t) = 
exp[— f(t)], a simple calculation gives 

g'(t) = -alog(l + pg{t)) +1 > 0, 

9{t) dy (2.23) 



*• Jo 



1 - alog(l +/jy) 



t. 



The first equation of (|2.23j) yields 1 + fJ,g(t) < e», whence sup t>0 g(t) < oo. 
On the other hand, the second equation of (I2.23P implies, for any fixed t, 
that g(t) (resp. f'(t)) is decreasing (resp. increasing) with respect to a. 
Thus the maximal and minimal solutions of (|2.22p take place respectively 
for a = p and a = £, and we have 

inf/'(t) ^expfr 1 ], VaG[£,p]. 
Consequently, by (|2.21|) . lim\e'(t) > 0, which contradicts null-recurrence, 

t— >oo 

and shows finally that if p > e -1 then necessarily £ < 1 (the announced 
transience). 

(ii) Consider now the case £ < 1. Letting s — > in (I2,14|) yields immediately 
Then rewriting (I2.14p in the form 

S£ * {s) _ (i-e)\p*(o)-/3*( s )] 

A s + s/3*(s) 

letting again s — > in (|2.24p and using l'Hopital's rule, we obtain 

/•oo 

e = \(\-tf I t(3(t)dt < oo. 



(2.24) 



The exact computation of £ proves to be a difficult project. Actually, since 
one can hardly expect more than approximate formulas, we shall present 
various results, both formal and concrete, some of them yielding bounds for 
£. 
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2.2.1 Formal approach 



Using the definition (|2.13p . it appears that the right-hand side member of 
(|2.14p can be analytically continued to the region 3?(s) < — Thus an 
analysis of singularities becomes theoretically possible, which should hope- 
fully allow to compute I. We roughly sketch the method, without giving an 
exhaustive presentation of all technicalities. 



Owing to the inversion formula ([2 
form 



we can rewrite (|2.14p in the functional 



1 

2in 



e st l3*{s)ds 



/iexp 



—A 
2iK 



e st ds 



;i +/?*(*)) 



R(<r) > 0. 



(2.25) 

Arguing by analytic continuation in (|2,25p , it is possible to prove that (5* (s) 
is a meromorphic function with real negative poles. Hence, f3(t) can be 
represented by the Dirichlet series 



P(t) = Cexp 



-At 



l + /3*(0) 



£ 

i>0 



me 



-ait 



(2.26) 



where C is a constant, the a^s form a sequence of positive increasing numbers 
satisfying 

ai> rrkoy Vi - ' 

and the iij's are ad hoc residues. In the ergodic case /3*(0) = oo and the first 
term in (|2,26p reduces to the constant C. Then, e(t) could be obtained by 
formal inversion of f3*(s). Alas, the computation becomes formidable and 
we did not get an exact tractable form (if any at all !) for £, since this is 
equivalent to compute Ui,ai,i > 0. 



2.2.2 Bounds and tail distribution 

Beforehand, it is worth quoting some simple facts. First, the value of I does 
solely depend on p, as can be seen by scaling in system (|2.1H2.2p . with the 
new functions 



Secondly, combining (I2.ip and (|2.2p leads to the inequality 

P(t)<n-\p(t), 
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which yields 



t < mini 1 



V/5 < OO. 



(2.27) 



In Section [2.11 we also could have considered the scheme 



70 (t) 
7k (*) 

?*(0) 



pe~ xt , 
dq k {t) 



+ 



dt 
pexp 

0, VA: > 



f > 0, 
t 

7& (* - y)dgk(y), 



(2.28) 



1 - 9ft 



which differs from (|2.10p only by its first equation, but this is a crucial differ- 
ence, corresponding in some sense to a fictitious function q~i(t) = 0, Vi > 0. 
Actually, this scheme produces a sequence of trees {L k ,k > 0}, with the 
property that the leaves of L k at level k never die. Its basic properties are 
the following: 

• the q^s form an increasing sequence of defective distributions; 

• for all k > 0, the tail distribution of q k dominates a defective exponen- 
tial distribution with density of the form akb k e~ bkt . Moreover, under 
condition (12. 51). we have 



lim cifc 

k— >oo 



lim 6^ = — 



and Qfc converges in L\ to the proper distribution p. 

The iterative scheme (|2.28p is convergent for all p, but the distributions 
Qk(t)j k > 0, are defective, their limit being proper if and only if p < e _1 . 
When p > e _1 , the limiting function p(t) remains defective and 



lim p(t) = lim lim qk(t) 



< 1. 



We shall derive bounds on £, in showing by induction that qk(t), for t suffi- 
ciently large, dominates an exponential distribution. The idea of proof will 
appear from the very first step k = 1. Actually, we have 

-e t\ 



90 (*) 



1 



A + /x 
7i (t) = /iexp 



-A(l 



|t + ^(e-*"-l) 
#o 
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and the Laplace transform 7*(s) has an explicit form, based on the formula 
(which involves the incomplete gamma function, see e.g. [9]) 



I(x,y) 



exp 



-xt + ye 



n=0 



n\(x + n) 



By scaling, for any constant c > 0, we have 



l^/x 



-,y 



exp 



o 



-xt + ye 



-ct 



9fc(x) > 0. (2.29) 



dt, (2.30) 



so that 



7i(«) = Y e * p 



s + X(l-£ ) X£o 
9q ' #o 



, R(« + A(l-£ )) > 0. 



The series in equation (|2.29p shows that 7*(s) can be analytically continued 
as a meromorphic function of s, with simple poles s n = —X(l—£o)—n,n > 0. 

Similarly, one checks easily the roots in s of 7*(s) + 1 = are simple, real 
and negative. Denoting them by —z n ,n > 0, we have the following 
Lemma 2.2. 

9l(s) - S (l + 7^)) 
is a meromorphic function of s, with poles at 0, — zq, — Z\, . . ., where 

A(l - £q) + n8 <z n < A(l - £„) + (n + 1)0 O , n > 0, 

wi/t £/ie more precise bounds 



Hence 



exp 



-A4 



< — A(l — ^o) — mm 



^o,/xexp 



-A4 



(2.31) 



(2.32) 



n>0 



where the residue r n of q\(s)e st at the pole z n ,n > 0, is positive and given 
by the linear relation 

TnzJ^r +1 = 0, 

CtS \s= — z n 



and 



Moreover, \2. 32j) yields 



n>0 



7i* (0) 
+ 7?(0)' 



h(l-e- zot ) <qi(t) <ii. 



(2.33) 
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Proof. Only the first part of (|2.3ip needs some explanation. It is obtained 
by checking that, for all y > 0, the first negative root in x of the equation 



px exp(— y)I(x, y) + Ox = 0, y > 0, 



satisfies 



xOq > — /xe y 



x *+(l + ^) +il e -j/ < o 



We shall prove by induction that p(t) dominates a reasonable exponential 
distribution. To this end, assume 



= (1 



which is in particular true for k = 0,1, as shown in lemma [2721 Then the 
calculus which led to fl2.33|> yields also 

4+i(l-e^ fc+1 ') <5fc+i(t), 



where (4+1, can be derived from (Ik, 6k) by the formulas 



Oik 



M cxp /-A4V 

Oik 



A(l-4) A4 



(2.34) 



1 + a k 

and is the first positive root of the equation 



• exp 



-A4Y 



-e k+x + A(i -4) A4 



6k 



+ 1 = 0. 



(2.35) 



Replacing an d zq in (I2.3ip by 9 k and 6 k +i respectively, one can prove the 
existence of 

lim (4, 6 k ) = (4, 6), 9 < oo, 4 < £ < 1, 

fe— >oo 

where < 9 < oo when £ < 1. 

The numerical computation of 4 is freakish in the ergodicity region (where 
the determination of the #fc's is a source of numerical instability), but proves 
very satisfactory for p^> e _1 . 

Next, instead of providing as in lemma 12.21 a stochastic ordering for all t, 
we get a tail-ordering, which has the advantage of achieving the exact value 
£ = 1 for p < e -1 . 
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Lemma 2.3. 



Qk(t) > o fc (l - e~ bkt ) + o(e~ M ), Vfc>0, (2.36) 
where the sequence satisfies the recursive scheme 



(2.37) 



(2.38) 



(-Xa k \ 

ak+ih+i = ^exp^— — J, 

b k+ i(l - ofe+i) = A(l - o fe ), 
wi/i ao = /V(A + /i) and bo = A + [i. 

Setting a — lim a& and 6 = lim in l\2.37\) . one has the limits 

fc— >oo k— >oo 

\ 

a = 1, b = — , if p < e~ , 
r 

a = x, 6 = A, if /0>e -1 , 
where x < 1 is the root of the equation 

xe x = -. (2.39) 
P 

In the course of the proof of lemma I2.3[ we will have to characterize positive 

f* (s) 

measures when then are defined from a Laplace transform of the form ^jWgj > 
as for instance in (I2.1ip . which a priori does not correspond to a completely 
monotone function, according to the classical definition of [6|. The following 
lemma does address this question and might be of intrinsic interest. 
Lemma 2.4. Let Q be a measure concentrated on [0,oo[, and its corre- 
sponding Laplace transform Q*(s) = f °° e~ st dQ(t), for any complex s with 
Sft(s) > 0. Define 



CO 



Ms) = f / pe-^ +st Ut, 
Jo 



oo(s) 



l + ^(s) 



Then co(s + fi) is the Laplace transform of a positive measure Aq [0, oo[. In 
addition Aq is a decreasing functional of Q, in the sense that, for all R> Q, 
R being Q- continuous, 

Afi < Aq. 
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Proof. For any complex number s with 5R(s) > 0, ip(s) = ip(s + //) can be 
viewed as the the Laplace transform of a positive measure U having the 
density Aie -(AQ(*)+^). Thus 

P{U <t}= f fie- {XQ{t)+ ^dt < 1 - < 1, (2.40) 
J o 

and the following expansion holds 

oo 

u( S + v) = Y,(-lN k+1 (s), (2.41) 

where ifj k (s) stands for the transform of the A;- fold convolution of U defined 
in (|2.40p . A function being uniquely determined — up to values in a set of 
measure zero — by the values of its Laplace transform in the region 3i(s) > 
[A, the first part of the lemma is proved. As for the monotony, one can 
differentiate the inverse of ([2.411) term by term (with respect to Q): since 
each term is multiplied by (— X) k , the resulting series is negative and the 
conclusion follows. ■ 

Proof of lemma [273[ Most of the ingredients reside in the integral represen- 
tation of 7 x (s) by means of formula (|2.29p . and we shall present the main 
lines of argument. 



Fix a number D, 6q < D < oo. Then (|2.30p yields the inequality 



7i(«) > ^i( s ) = £ exp (~^ ) 2 



s + A(l — do) ^ a o 



(2.42) 



whence 



Setting 



SQi is) = — — > , . , , Ms > 0. 



n(g) = ^exp(-^j A(i - TT - , (2.43) 



— Aao^ 1 

ao) + s " 

and isolating the first term is the power series expansion of (|2.42p by means 
of (|2.29p . we obtain after a routine algebra 



V>i(s) _ u(s) | w(s,D) 



l + ^i(s) u(s) + l D 
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Remarking the first pole of w(s, D) is the root of u(s) + 1 = 0, we can use 
lemma 12.41 and a term by term inversion of (|2.44p to obtain 



qi(t)>a 1 (l-e- b ^ t ) + ^^ 



D 1 

where the couples (ai,&i), (ao,&o) satisfy system (|2.37p . 

At step k = 2, one would introduce a constant, say D\, and repeat the same 
procedure to obtain (|2.36p . It might be useful to note that it is not possible 
to take D = oo, since this would create an atom at t = 0, in which case 
lemma [231 does not work in general. 

At last, it is a simple exercise (therefore omitted) to verify the existence 
of (a, b) = limfc_ >00 (afc, 6^), given by (|2.38p . The proof of the lemma is 
terminated. ■ 

Now, to conclude the proof of theorem 12.11 it merely suffices to note that, 
by 4139]), 

P' 1 - P~ 2 < x < £ < p- 1 , V/j^e -1 . 



Subsidiary comments The method of schemes to analyze nonlinear oper- 
ators in a probabilistic context is extremely powerful (see e.g. [3] for problems 
related to systems in thermodynamical limit), and in some sense deeply re- 
lated to the construction of Lyapounov functions. Up to sharp technicalities, 
the schemes (I2.10P and (12.28H can be exploited to derive precise information 
about the speed of convergence as t — > oo, for any p, < p < oo, and 
when pushing exact computations slightly farther, one perceives underlying 
relationships with intricate continued fractions. Finally, we note that the 
question of transience could be studied from a large deviation point of view, 
by considering e(t) as the member of a family indexed by the parameter 
(p — e _1 ) — see in this respect section [4) 



3 Some stationary distributions 

In this section, we derive the stationary laws of some performance measures 
of interest when the system is ergodic, i.e. p < e _1 . Incidentally, note that 
the only process studied so far, that is the number X of vertices attached 
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to the root, behaves like the number of customers in a m/g/oo queue with 
arrival rate A and service time distribution p, so that 



lim P{X(t) = k} = e~ r1 —, \fk > 0. 

Another point worth mentioning is that, as for the model in |16| . the Markov 
process G(t) is reversible and hence has an explicit invariant measure. To see 
this, notice that at each vertex v, leaves are added at rate A, and removed 
at rate fJ-rj(v), where rj(v) stands for the number of leaves attached to v. 
Therefore, the stationary probability of some configuration G is 

n N G 

where K is a normalization constant. The ergodicity of the Markov process 
G(t) is then equivalent to the convergence of the series X^ 7r (^')> where the 
sum is taken over all admissible trees G. However, while counting Catalan 
trees as in |16j is not that difficult, the combinatorics is more involved in the 
present setting, and this direction will be pursued no further. 



3.1 Volume of the tree 

Let N = lim^oo N(t), where N(t) introduced in section Q] stands for the 
volume of G(t) and the limit is taken in distribution. 

Theorem 3.1. When p < e _1 , the distribution of the stationary volume N 
is given by 

P{N = k} = -^-p\ (3.1) 
r k\ 

where r is given by h2. 6\) . Moreover, the mean value of N is given by 

EN = — !— . (3.2) 
1 — r 

Proof. We proceed as in lemma [2TT| saying that the number of vertices in the 
tree at time t is equal to 1 plus the numbers of vertices in all the descendants 
that have appeared in [0, t] and are not yet dead. The construction mimics 
the former proposed for the process X v (t): the volume of a subtree rooted 
at some vertex v is distributed as N(t) for t < r v . 
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For any complex number z, \z\ < 1, we have therefore 



fc=0 



e- xt (\tf 



dx , 



y N(x)t {x < T} 



exJ\E[j\z N W - l)l {x < T} dx] \, (3.3) 



and, letting t — > oo, 



Ez N = zexpi AE 



(z N ^-l)dx 



pz 



exp< AE 



z N ^dx 



L J0 



• (3-4) 



It is easy to write a renewal equation similar to (|2,2p . namely 

Ez N ^ = E[z N (H {t<T} ]+ f Ez N ^dp(y), 

Jo 

which, after setting 4>(z,t) = Ez N ^ and 4>(z,t) = E[z N ^1\-j t < T -i], and taking 
Laplace transforms with respect to t, yields the equation 

0*0, s) = 4>*{z,s)+4>*(z,s)sp*(s). 

Then, as in the case of (|2.8[) . using the boundedness of z N ^\ we get 

i r°° ~, , . 1 



(z) = lim <j)(z, t) 

t— yoo 



<fi(z,t)dt = — E 
m ./ m 



z N ^dx 



Then (|3.4p can be rewritten as 

r<f>(z) = pz exp{rcj)(z)\ , 



(3.5) 



and hence rcj)(z) = C(pz), where C stands for the classical Cayley tree 
generating function (see e.g. |17j). Using the well-known series expansion for 
C (which follows from Lagrange's inversion formula), we get (|3.ip . since 



1 ^ k k - x 



r < k\ 



(pz) k . 



The mean (|3.2p is obtained by differentiating (|3.5p with respect to z and 
taking z = 1. ■ 
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The analysis of the asymptotics of (|3.ip with respect to k confirms an in- 
teresting change of behavior when p = e~ l . Indeed, for k sufficiently large, 
Stirling's formula yields 

pw = k}= l - k —o k « 1 * = ± 

r k\ r^kk k e~ k r^ k l 



Moreover, a straightforward Taylor expansion of (|2.6p gives the following 
estimate of EN, as p — > e~ 1 : 

EN 1 



1 - r ^2(1 - P?) 

Thus, while all moments of N exist for p < e _1 , there is no finite mean 
as soon as p = e _1 . We note in passing that this phenomenon appears 
sometimes in branching processes and can be viewed as a phase transition 
inside the parameter region, as already remarked in [5j. 

3.2 Height of the tree 

Let H = lim.t_j.oo H(t), where H(t) introduced in section Q] stands for the 
height of G(t). The distribution of H is given by the following theorem. 
Theorem 3.2. 

1. For p < e~ 1 , the following relations hold: 

P{H = 0} = e~ r , (3.6) 
P{H>h + l} = 1 - exp[-rP{F > h}], V/i>0. (3.7) 

2. If p < e _1 , then there exists a positive constant 9{r, 1), such that 

2h 

P{H>h} = 9{r,l)r h+1 + o( K ^—^j (3.8) 

where the function 9{r,x) is the locally analytic w.r.t. x solution of the 
functional Schroder equation 

9{r,l-e- rx )=re{r,x), 

subject to the boundary condition 

f)9 

S M) = 1. (3.9) 
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3. Whenp = e~ x , 



P{H >h} = - + O 



log h 



(3.10) 



Proof. As in the previous proof, one writes the height of the tree at time t 
is less than h + 1 if, and only if, all the descendants that have appeared in 
[0, t] are either dead or have a height smaller than h: 

P{H(t)<h + l} = ±^^{J^[l-P{H( X )>h,x<r}]} k 
= exp 

{-^J P{H(x) > h,x< r}cte}. 

Letting t — > oo and arguing as in theorem 13. 1| we can write 
P{H <h + l} = exp[-rP{# > h} , 



which proves (|3.7p . On the other hand, equation (|3,6p is immediate, since it 
is in fact a plain rewriting of (|2,3p . 



To prove the remainder of the theorem, let do = x, where x is a positive real 
number, and consider the sequence 



d. 



h+l 



1 



-rd h 



, h = 0,l,... 



(3.11) 



When x = 1, note that we have exactly dh+i = P{H > h}. The question 
that faces us now is to compute and to estimate the iterates of an ana- 
lytic function, in the circumstances 1 — e~ rx . This subject concerns a wide 
branch of mathematics (including functional equations, automorphic func- 
tions, boundary value problems), and it has received considerable attention 
since the nineteen twenties. We shall employ classical arguments without 
further comment, referring the interested reader to e.g. [12j and [lj for a 
more extensive treatment. 

A Taylor expansion up to second order in (|3.1ip gives 

rd h > 4+i = 1 - e~ rdh >rd h - r -^, (3.12) 

which implies that r~ h dh is a decreasing sequence with lim \. d^ = (that 

h— >oo 

we already knew!) and 

d h < xr h . (3.13) 
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As h — > oo, the asymptotic behavior of dh has a twofold nature, depending 
on whether r = 1 or r < 1. 



Case r = 1. This is the easy part. Writing 



' ' ' +\ + 0{d h ) 



d h+1 l-e- d h d h 2 
we get immediately dh = O ( jr J . and hence 



+ 0(5> fe ) =|+0(log/ i ) ! 



4 2 

fc=0 



which leads to (13.10|) . 

Case r < 1. The analysis is less direct. From (I3.12|) and (13. 13|) . we infer 
that, when h — > oo, r~ h dh has a limit denoted by 9(r,x), with 



< r~ h d h - 9(r,x) < 



9 h 

rx r 



2 1-r 

First let us show that 9(r,x) is strictly positive. Indeed, 

d ^=x\{{l-^ m {r,x)), (3.14) 

where, Vm > 0, the quantity ip m (r,x) = 0{r rn+l ) is an analytic function of 
the pair of real variables (r, x) in the region [0, l[x[0, A], with < A < oo. 
Hence, as h — > oo, the infinite product in (13. 14ft converges uniformly to a 
strictly positive value, Vx > 0, so that 9(r,x) is also analytic of (r, x) in the 
aforementioned region. To summarize, 

lim r~ h d h = 9(r,x) > 0. 

h— >OD 

The pleasant fact is that 9, taken as a function of x, satisfies the so-called 
Schroder equation 

9{r,l-e~ rx ) =r9(r,x). (3.15) 
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While it is clear that 0(r, 0) = 0, ([3.15P does not impose any constraint on 
d9 

— (r, 0). However, it is easy to show by induction that d) l (x) (where the 
dependency on x is for a while explicitly written) satisfies 

^ =r\Vh>0, 

OX \x=0 

and thus condition (|3.9p also holds for 9. To conclude the proof of f|3 . 8[) . it 
suffices to choose x = 1. ■ 



Remark We have taken the variable x on the positive real half-line to 
get sharper bounds, e.g. ()3. 13j) . Actually, arguing as above, it is immediate 
to check that 6 has an analytic continuation in the complex x-plane in a 
a neighborhood of the origin. In this respect, without going into a full 
discussion, we mention the relationships with automorphic functions and 
boundary value problems, which would allow integral representations. For 
our purpose, simply writing 



0(r,x)=J29ix\ 0o = O,0i = l, 



i>0 

we see that all the 0j's can be computed recursively. Furthermore the itera- 
tion of (I3.15P yields 

0(r,d h )=r h 9(r,x). 

from which we obtain 

d h = uj(r,r h 6{r,x)), 

where u(r, x) denotes the inverse function of 6 with respect to the variable x 
and satisfies the functional relation 

1 — exp{— ru(r, y)} = u)(r, ry). (3.16) 

We have 



w ( r >y) = ^2^iy\ wo = o, wi = i, 

and again the Wj's are obtained recursively. 
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4 Scaling and limit laws in the transient case 



In this section, we present some limit laws for N(t) and H(t), which are 
especially of interest when the system is transient. Beforehand, for every 
integer k and all t > 0, we define the quantities 

'X k {t) d = f #{v G G{t) : h(v) = k}, 

< oo 

Y k (t) d = f J2 X ^{t<r}- 

j=k 

So, Xk{t) stands for the number of vertices at level k in the whole tree at 
time t. 



4.1 Scaling for N(t) in the pure birth case fj, — 

Lemma 4.1. When fi = 0, EX n (t) has the explicit form 

(At)" 



EX n (t) = 



(4.1) 



Proof. Since 

P{X n (t + dt) = X n {t) + l\X n {t), X n _i(t)} = XXn.^dt + o(dt), 



we obtain 



-EX n (t) = XEX n ^(t), n>l, 



EX (t) = l, 

and the result is immediate by induction. ■ 
Theorem 4.1. When \i = 0, the expected volume at time t is given by 

EN{t) = e xt , (4.2) 

and 

Arm 

(4.3) 



lira J^2r =Exp(l), 



t^oo EiV(t) 

where the limit is taken in distribution and Exp(l) denotes an exponentially 
distributed variable with parameter 1 . 
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Proof. Equation (|4.2p is a mere consequence of lemma 14. 1| since 

oo 

EN(t) = Y J EXn(t) = e xt . 

n=0 

Let now (f>(z,t) = Ez N ^ , for z complex with \z\ < 1. We start from equa- 
tion (|3,3p . in which we take r = oo. Then the following relation holds: 

<fi(z, t) = z exp|A J (<p(z, x) — l)dx 

Differentiating with respect to t yields 
d 



-<j)(z,t) = \[4>{z,t)-l]0(z,t), 



whence 



<K*,t) ' 



where K does not depend on t. Since (p(z, 0) = z, we deduce K = z 1 — 1, 
and finally 

4>{z,t) 



l 



1 + [z- 1 - 11 



.At ' 



The Laplace transform of e xt N(t) is, for K(s) > 0, 

Eexp{-se~ A *7V(t)} = (/>(exp{-se~ At }, i 



1+ [exp{se- Ai } - l]e At ' 



so that, letting t — > oo, 



Urn Eexp{-se~ A *7V(t)} = - 



Now (|4.3p follows directly from Feller's continuity theorem (see [6]). 



4.2 An ergodic theorem for H(t) 

The key result of this section concerns the height of the tree and is formulated 
in the next theorem. 



Let 



1 / x def S 

Ms, c) = - + log 

c 



\(l-sp*(s)) 



, »(s) > 0. 



(4.4) 
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Theorem 4.2. With probability 1, 

Urn*®.* 

where 5 >0 is uniquely defined from the system of equations 
In f/ie ergodic case, 5 = 0. 

The proof is constructed around the three forthcoming lemmas. 
Lemma 4.2. Define the events 

Hit) 1 „ r, £T(t) 



A c = (liminf > c), 5 C = (limsup < c\. 



T/ien P{A C } = or 1 and P{i? c } = or 1. /n oi/ier words, A c and B c satisfy 
a zero-one law and can only be trivial events (i.e. sure or impossible) . 

Proof. Fixing an arbitrary to, with G(to) = Gq, we want to show that A c 
does not depend on Go- For this purpose, consider the random process 
G'(t) C G(t) constructed as follows: for t < to, it consists only of the root, 
and for t > to it contains exactly that part of G grown from the root after 
time tQ. Then the probability that A c holds for G' is clearly equal to the 
probability that A c holds for G without conditioning. In other words, since 
Ha{t) > we have 

P{A C | G(t ) = G } > P{A C }. 

Then basic properties of the conditional expectation yield 

E[P{^ C | G(t )}] = P{A C }, 

so that 

P{A C | G(t ) = G } = P{A C } (4.5) 

for any Go- On conditioning with respect to G(to), G(t\), . . . , G(tk), for any 
arbitrary increasing sequence of times t^, we see (|4.5p still holds. Hence, the 
assertion for A c is a direct consequence of the zero-one law for martingales 
(see e.g. [ID]). 



Quite similarly, if the event B c holds for G, then it is also in force for any 
subtree rooted at a vertex of Go, which reads 

P{B C | G(t ) = G } < P{B C }. 

The lemma is proved. ■ 
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Lemma 4.3. 

(i) If, for some integer n and real number c > 0, E[Y n (n/c)] > 1, then 

P{A C } = 1. 

oo 

(ii) If for some n and real number c > 0, E[Xi cn (kn/c)] < oo, then 

k=0 

P{B C } = 1. 

Proof. 

For the sake of brevity, let J n (k) denote the time interval [kn/c, (k + l)n/c]. 

(i) Consider a standard branching process (,k,k>0, endowed with an off- 
spring distribution equal to that of Y n {n/c). From the condition in 4.3(i), 
this process has a probability of non extinction which is strictly positive and 
will be denoted by y(n, c). The key point is that can be viewed as defining 
a subtree G' C G such that HQ/(kn/c) > kn. Then P{A C } > y(n,c) > 0, 
and we have P{^4 C } = 1 by lemma [4~2l 

Indeed, to build such a subtree G' , we associate with each generation of 
a set of vertices Sk C G(kr), such that = \sk\- 

Let So = {^o} and, for each v € Sk,k > 0, let G v be a subtree rooted at v 
and born during J n (k) (by convention G v is empty if v dies). We put 

Sk+i = [J {v £ G v : d(v',v) > n}. 

v£s k 

This construction produces the desired tree, since the volume of each set 
belonging to the above union is exactly distributed as Y n (n/c), and because 
Sk consists of vertices located at a distance at least kn from the root. 

(ii) Let ak = o{k) be a sequence of non-decreasing positive integers. Then 
for any fixed integer ko, we have the inequality 

P(B C ) > P{ sup H(t) < (k + l)n + a k ,Vk > k \. 

We observe the height of the tree decreases at a rate not faster than fj,, so 
that, given the event {H((k + l)n/c) < (k + l)n}, the supremum of H (.) on 
the interval J n (k) is bounded by 

sup H (i) < (k + l)n + 7r(/m/c), 
ieJ„(fc) 
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where ir(x) denotes a Poisson random variable with rate x. Thus we have 

P(-Bc) > P{H(kn/c) < kn,Vk > k } f[ P{n(fJ,n/c) < a k }, (4.6) 

fc>fco 



and we will show that the right-hand side of (|4.6j) can be rendered positive. 
First, we remark that 

P{H(kn/c) <kn,Vk>k } > 1 - ^ P{H(kn/c) > kn} 

k>ko 

> 1 — E[X kn (kn/c)} — > 1, as ko — > oo. 

fc>feo 

Secondly, we choose the sequence 

a k = j, Vk G - l)/2 + 1, j(j + l)/2],Vj > 1, 
which consists of blocks of repeated integers satisfying a k = O^k 1 / 2 ). 

Setting v = fin/c, the product in (|4.6p will be positive, provided that the 
following sum is finite 

£pw„)>« fc } < E^^f 1 --)" 

a k- ^ &k' 

k>ko k>ko 

V 3 ( \ ~ 1 / 1/ N — 1 



< y e "V — tttI 1 --) < ^( 1 -^) <°°> 
- (j-l)!V j7 - V 



and hence (i^ follows from the zero-one property of B c . 

The proof of the lemma is concluded. ■ 

Lemma 4.4. For any 3?(s) > 0, let (p k (s) and <pk( s ) be the Laplace trans- 
forms 

POO 

<p k (s) = / EX k (t)e- st dt, 
Jo 

poo 

<p k (s) = f / E[X k (t)l {t < T} ]e- st dt. 
Jo 



Then 



X k (l-sp*(s)) 
Ms) = 3Ti 



<Pk( s ) 



fc+1 



A fc (l - sp*(s)) kA 
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Proof. It is not difficult to check the following relations 



EX fc (i) = I E[X k ^{y)t {y < T} ]Xdy, k>l 
Jo 

EX k (t) = E[X k (t)l {t < T} ] + / EX k (t-y)dp(y), 

Jo 

with the initial condition EAo(i) = 1. Actually, the first equation follows 
from an argument already employed before. Namely, the number of vertices 
at level k are the direct descendants of vertices at level {k — 1) still alive at 
time t, remarking that each such descendant on [0, t] appears independently 
at rate A. The second equation is a straight regeneration relation. Therefore, 

<fk(s) = k > 1, 

s 

ip k (s) = (f k (s) + ifk(s)sp*(s), 
whence, since (fo(s) = 1/s, 

A(l - S p*( S ))^ fc _i(s) _ A fc (l - s P *(s)) k 



s k+i 



and the result follows. ■ 
We are now in a position to prove theorem 14.21 

Proof of theorem The proof is split into two parts, each one correspond- 
ing respectively to criteria (i) and (ii) of lemma [4~3l 

First we shall find the largest c, denoted by Ci n f, ensuring criterion (i) of 
lemma 14.31 is fulfilled. Applying the results of lemma 14.41 and the inversion 
formula (12.71). we have 



i ro+ioo 

E[Y n (n/c)} = ^ — / ^{sy^ds 



j=n 
1 

2iK 



a+ioo 



ct— too 



\(l-sp*(s)) /c 



s 



l-sp*(s) 
s — A(l — sp*(s)) 



1 r +l °° e nb ^[l-sp*(s)] ds 



ZinJa-ioc s - A(l - sp*(s)) 
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intheregionW = {u > 0, a > A(l— ap*(a))}, where b(s,c) has been defined 
in (I4.4p , When the system is ergodic, it is immediate to check the region IA 
coincides with the complex half-plane !R(s) > 0. On the other hand, in the 
transient case, the equation 

s = A(l- sp*{s)) 

has exactly one root, which is real and belongs to the open interval ]0, A[. 
Computing the residues of the integral (|4.7p (by shifting the line of integra- 
tion to the left, after analytic continuation of p*(s) to the region a = — e, 
for some e > 0) is a tedious task, in particular due to the pole of order n 
at s = 0. We will rather proceed by a kind of saddle-point approach (see 
e-g. 0). 



The form of the integrand in (14, 7ft shows that, as n — > oo, the boundedness of 
E[Y„(n/c)], depends primarily on the value of the modulus of b(s, c). In fact 
one can see precisely that E[Y n (n/c)], for each fixed c, does not tend to zero 
iff the minimum of b(s, c) is non-negative at any possible real saddle-point 
s 6W, where 

^M = 0j SGU. 

OS 

It follows that Cj n f is the unique real solution of the system 

b(s,c mi ) = db{s : Cmi) =0, seU. (4.8) 

OS 

Without presenting a detailed discussion, we shall simply stress that in the 
real plane (s, y) the curves 



y = s/c and y = — log 



\(l-sp*(s)) 



are tangent (resp. intersecting, non-intersecting) for c = c; n f (resp. c > Ci n f, 
c < c inf ). 

As for the second part of the theorem, the question is to find the value Cs Up , 
equal to the smallest positive number c satisfying criterion (ii) of lemma |4T3| 
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which implies the finiteness of the quantity 



fc=o fc=o • y < T - JO ° 

A(l-sp*(s)) /c 



1 r +io ° (is r /W1 ^ n ^" 1 

-e° 



2i7T ./ a—ico S 



i /-cr+ioo J 



where (|4.9p holds in the region V = {cr > 0, cr > A(l — crp*(cr))e <r / c }. 

Clearly, the existence of the last integral in (|4.9|) . as n — >■ oo, amounts again 
to find the sign of 3?(6(s,c)), for s £ V, Arguing exactly as above, one can 
find at once the equality 

— f x 

^sup — Qnf — 0. 

When the system is ergodic, lim^o b(s, c) = log Am = logr < 0, which 
yields 8 = as might be expected. 

The proof of the theorem is concluded. ■ 

As a by-product, we state the following corollary, of which the almost sure 
convergence part has been derived in [4j [15] through different and less terse 
methods. 

Corollary 4.1. In the pure birth case fi = 0, almost surely and in L\, 

li m E^L = Ae. (4.10) 

t-S-oc t 

Proof. Instantiating equation ()4.1j) in criteria (i) and (ii) of lemma l4.3l yields 
directly the first part of (|4.10|) . On the other hand, it is immediate to check 
that the function EH(t) is superadditive (this would be not true for [i ^ 0), 
namely 

EH{s + t) > EH(s) + EH(t), 
so that, by a variant of Kingman's theorem (see |10j). the limit lim — 

t— >oo t 

does exist. Then the convergence in L\ will follow if one can show 

EH(t) < At, Vt > 0, 
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for some positive constant A. Using the obvious inequality 
P{H(t) >k} = P{X k (t) > 0} < EX k (t), 

we can write 

oo fe oo f\ + \k 

E/r (t) = £ >*}<£i + £ ^r- 

fc=l k=l k=k +l 

Then, taking ko = \Xet] and using Stirling's formula, we obtain 

rrr . , , {Xt) k ° +1 Xt \k (Xt) Xet+1 e 

EH(t) < ko + i—t tt) <Xet+ — v ; ; 

W " {k + l)\f-Ak + lJ ~ (Xt)^V2^XYte-l 



< Xet + 



k=0 



2vr(e - 1) 



5 Extension to the multiclass case 

The extension of the results of section [2] to models encompassing several 
classes of vertices is very tempting, although not quite evident. We solve 
hereafter a case where the birth and death parameters depend on classes in 
a reasonably general way. 

Let C be a finite set of classes. Then the multiclass Markov chain Gc has 
the following evolution. 

• At any given vertex of class c, a new edge of class d G C can be added 
at the epochs of a Poisson process with parameter A cc / > 0. 

• Any leaf attached to an edge of class d and having an ancestor of class 
c can be deleted at rate fi cc > > 0. 

• The root vq of the tree belongs to class c G C, say with probability tt c , 
with X^cec 11 c = 1> a lbeit these probabilities will not really matter in 
our analysis. 

Let p cc i be the lifetime distribution of a vertex of class d which descend from 
a vertex of class c. Also, X c (t) will denote the distribution of the number 
of direct descendants of a vertex of class c. The following lemma is the 
analogous of lemma 12.11 
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Lemma 5.1. The lifetime distributions p cc > , c,c' G C satisfy the following 
set of equations. 

P{X c (t) = 0} = exp{-^A cc , I (l- Pcc ,(x))dx}, (5.1) 
P{X c (t) = 0} = J_^l + /"* P{Xc (t -y) = 0}d Pbc (y), V6 G C, 

Pbc Clt Jo 

(5.2) 

with the initial conditions p cc '(0) = 0, Vc, c' G C. 

Proof. Details are omitted, as it suffices to mimic the proof of lemma 12.11 
Note however that the dependency with respect to d disappears surprisingly 
enough in (I5.2p . Indeed, while the lifetime of a vertex depends on the class of 
its direct ascendant, the distribution of the number of its descendants merely 
depends on its own class. ■ 

In the setting of this section, it is actually not easy to come up with a natural 
explicit extension of theorem 12.11 However, the following theorem provides 
a necessary and sufficient condition for ergodicity. 



def A cc / 



The following notation will be useful in the theorem: 

Vc, c G C 

c'&C 

We will also denote by p > the Perron- Frobenius eigenvalue (see [8]) of the 
non-negative matrix (/°cc') c c / gC - 
Theorem 5.1. 

1. The Markov chain Gq is ergodic if, and only if, the nonlinear system 

yc = ^2 pcd ex p{^}' Vc G C ' ( 5 - 3 ) 

dec 

has at least one real-valued (and obviously non-negative) solution. In 
this case, the mean lifetime of a vertex of class c with an ascendant of 
class b can be written as 

m bc = — exp{r c } , (5.4) 

P-bc 
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where the r c form the smallest solution of \5. 3\) . that is r c < y c , Vc G C. 
Note that h5.3\) implies that r c > p c . 

2. A simple sufficient condition for ergodicity is 

Pc<~, Vc G C (5.5) 

in which case r c < p c e. 

3. A simple necessary condition for ergodicity is 

p<~. (5.6) 

e 



Remark Before stating the proof of the theorem, it is worth pointing out 
that equation (I5.3P may in general have several real solutions (as in dimension 
1). Therefore, there is no guarantee that the solution y c is the correct one. 
However, its sole existence proves ergodicity and (I5.4j) . 



Proof. Assume first that Gq is ergodic. Then, as in theorem 12.11 we let 
t — > oo in (15.1I)-(15.2H to obtain the relation 



- — = exp{- S~] A cc /m cc /|, 



P-bcmbc ,_ ( , 



which in its turn yields (|5,3p and (|5.4p . just choosing 



c'ec 



As for the proof of sufficiency in item [H we introduce the following modified 
version of scheme (|2. 10[) : 

Pcc'-fi{t) =1, t > 0, 

a c;k+1 (t) = exp|-V]A cc / / (1 - Pcc';k{y)) d y\' k >°i 

1 cec Jo t J (5.7) 

a c -,k{t) = — dpbc ^ + j a c . k {t-y)dp bc . k {y), k > 1, 

P'bc 

[Pcc';k(0) = 0, k>l. 
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Then, for any b, c G C, we have 

* c; i(t) = 1 and Pbc;1 (t) = 1 - e - "** < p bc - (t). 



a, 



Here again, the positive sequences {p C c';k(t)'> k > 0} and {a C ;k(t)] k > 0} are 
uniformly bounded and non-increasing, for each fixed t > 0. Consequently, 

p cc /(t) = lim\p cc /. fe (t) and a c (t) = lim\ a c;fc (i) 

form the unique solution of (|5.ip - (|5.2j) . uniqueness resulting from the Lips- 
chitz character of equation f|5 . 2[) . 



Letting m cc t.^ denote the finite mean associated with each distribution p cc '-k 
and 

dcf . 
r c;k — / J ^cc' m cc';k > 

c'eC 

we can write the following recursive equation 

r C ;k+i = ^2 Pec' exp{r c / ;fc }, Vc G C. 
dec 

The Pcc'-k 's are decreasing sequences, and hence the r c; fc's are non-decreasing, 
with r c -Q = 0, Vc G C. If (j5.3j) has a solution, then the relation 

r c ±+i ~ Vc = ^2 Pcd [ ex P{ r c';fcl ~ ex P{yc}] 
c'eC 

yields r c -k < y C) for all c G C. Therefore, each sequence r c; fc converges as 
A; — > oo to a finite value r c < y c , and Gq is ergodic since, by (|5.4|) . the m cc /'s 
are also finite. When (I5.5P holds, the same line of argument shows that the 
sequences r c -k are non-decreasing and bounded from above by p c e. 

Finally, to prove (|5.6|) . we use the following inequality (see |8J), valid for any 

x c > 0, c G C: 

, Pcc' x c' 
p < max > . 

c'eC c 

When the r c 's satisfy (|5.3p . the choice x c = exp{r c } implies 

1 



p < max 

C6C 



^ p cc > exp{r c /} exp{-r c } 
c'eC 



= max 


r c exp{-r c } 


cec 





< 

e ' 



which concludes the proof of the theorem. 
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It is possible to extend the results of section [3] to the multiclass case. We will 
only sketch the proofs in what follows, since they are very similar to their 
single class counterparts. At time t > 0, let N c d(t) be the number of vertices 
of class d inside a tree, the root of which is of class c. Then, as in proof of 
theorem 13. 1| if z c is a complex number such that \z c \ < 1, Vc G C, 



[n 

dec 



,N cd (t) 



expjEMAlI^- 1 ) 11 

^c'ec Jo dec 



i x < T cc'} 



dx 



Assume the system is ergodic, let N C( i = lim^oo N c d(t) and 

*c(2) = E[lK- 



dec 



Then computations similar to the ones in theorem 13.11 yield 

(f> c (z) = z c expl ^2 Kc' , mcc'(4>c'(z) - 1) >, ceC. 
^c'ec ) 



(5.8) 



Unfortunately, no closed form solution is known for cf) c from this equation. 
It is however possible, as for (|3.2p . to write down a system of equations for 
the expectations of the N^s. 

E[N cd ] = l{ c =d} + Pec' exp{r c /}E[A r c / d ] . 

c'ec 



This system admits of a non-negative matrix solution if, and only if, the 
Perron-Frobenius eigenvalue of the matrix 



M = (/9 cc /exp{r c /}) 

V / c,c 



c'ec 



is smaller than 1. A simple necessary condition for this to hold is (15. 5ft . 

Finally, the same line of argument allows to extend (|3.7p . If H c is the height 
in stationary regime of a tree which root is of class c G C, then 



P{H c = 0} = e~ r % 
P{H C > h + 1} = 1 - exp 



^ p cc , exp{r c /}P{# c / > h} 

c'ec 



, Vh > . 
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